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Abstract 

In this article we investigate traveling wave solutions of a nonlinear differential equation 
describing the behaviour of one-dimensional viscoelastic medium with implicit constitu¬ 
tive relations. We focus on a subclass of such models known as the strain-limiting models 
introduced by Rajagopal. To describe the response of viscoelastic solids we assume a non¬ 
linear relationship among the linearized strain, the strain rate and the Cauchy stress. We 
then concentrate on traveling wave solutions that correspond to the heteroclinic connec¬ 
tions between the two constant states. We establish conditions for the existence of such 
solutions, and find those solutions, explicitly, implicitly or numerically, for various forms 
of the nonlinear constitutive relation. 

Keywords: traveling waves, viscoelasticity, strain-limiting model, implicit constitutive 
theory. 


1. Introduction 

The present paper is concerned with the dynamics of a viscoelastic medium investi¬ 
gating the traveling wave solutions of the equation 

T xx + vT xx t = g(T) t t, (IT) 

where T(x, t ) is the Cauchy stress at point x and time t , g(-) is a nonlinear function and 
v > 0 is a constant. Equation ED is a one-dimensional nonlinear differential equation in 
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T resulting from the equation of motion and a constitutive equation relating the stress, 
the linearized strain and the strain rate. 

As opposed to the classical models in mechanics, the strain can be written as a 
function of the stress, rather than expressing the stress in terms of the kinematical 
variables. This idea is due to Rajagopal |ll, 12]] > who introduced a generalization of 
the theory of elastic materials by suggesting implicit models allowing for approximations 
where the linearized strain is a nonlinear function of the stress. A series of pap ers on 
such implicit theories has been published recently (see e.g. fl, fl, fl, H, M)- The 
advantage of this new idea is that it allows for the gradient of the displacement to stay 
small so that one could treat the linearized strain, even for arbitrary large values of the 
stress. In this work we focus on four different such models, and we reconsider them in 
the context of viscoelasticity. We also look at models with quadratically and cubically 
nonlinear constitutive relations although they do not behave as expected for large values 
of the stress. 


There are numerous models introduced by Rajagopal in [ll| with implicit constitu¬ 
tive relations between the stress and the strain including models for elastic fluids, in¬ 
elastic materials and non-hyperelastic materials. Following these models, various forms 
of nonlinear constitutive relations have been studied in different contexts. For example, 


Kannan, Rajagopal and Saccomandi fl worked on the elastic case with a polynomial 


• fl, 


type nonlinearity (see Section [ 2 ] for more details). Bulicek et al. [lfl, on the other hand, 
considered the static case with a more general nonlinearity (see Section [ 2 ]) and presented 
the first existence result in a three-dimensional domain. 

For viscoelasticity, much less is done in the literature. As explained by Muliana, 
Rajagopal and Wineman in [jfl], force, and hence the stress, is the cause for deformation, 
hence for the strain. Because of this the strain should be described in terms of the stress 
or its history than vice versa. The motivation for this idea is that in the classical elasticity 
theory, there cannot be a nonlinear relationship between the linearized strain and the 
stress, which, in fact, is observed in some experiments (see e.g. [19)1 , [lfl ). The fracture of 
brittle elastic bodies is another possible application area for such implicit theories, where 
one can obtain bounded strain at the crack tip due to the possibility of having a nonlinear 
relationship between the linearized strain and the stress (see [jj]] for details). Muliana 
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l. H 


et al. m developed a quasi-linear viscoelastic model where the strain is expressed as an 


integral of a nonlinear measure of the stress. Rajagopal and Srinivasa in 


17| proposed a 


Gibbs-potential-based formulation for the response of viscoelastic materials in this new 


class. Also Rajagopal and Saccomandi m investigated viscoelastic response of solids, a 
one-dimensional version of which is the one we study in this work, namely 

yB + i^D = /3 0 I + /3 iT + /3 2 T 2 , (1.2) 

where 7 and v are nonnegative constants, fa = fa(h, h, h), (i = 0,1, 2), I\ = trT, / 2 = 
±trT 2 ,/ 3 = itrT 3 , B is the left Cauchy-Green stretch tensor and D is the symmetric 
part of the gradient of the velocity field. As they explain, this model includes as special 
subcases; models for a very general new class of elastic and viscoelastic bodies (e.g. 
Titanium and Gum metal alloys), as well as the Navier-Stokes fluid model (see e.g. [12| ). 
Linearizing the strain in this model reduces CE2D to 


e + vet — A)I + /3iT + /3 2 T 2 , 


(1.3) 


where e = |(Vm + Vit T ) is the linearized strain, and u(x,t ) is the displacement. 

We study m in one space dimension with a general nonlinear right-hand side (see 
m)- We are interested in analyzing the conditions on the nonlinearity g(T) when 
traveling wave solutions of the form T(£) with £ = x — ct, where c represents the wave 
propagation speed, exist for two constant equilibrium states at infinity. We find the 
solutions analytically (implicitly or explicitly) or numerically. More precisely, we will 
first look at the quadratic and the cubic cases for which we are able to solve the problem 
analytically and obtain explicit or implicit solutions. After that we will study four 
nonlinear models, namely Models A, B, C and D (see Section [2]), and we will either 
express the solution implicitly, or obtain it numerically if it is not possible to find an 
analytical solution. Our work seems to be the first such treatment in the literature of 
strain-limiting viscoelasticity. 

The propagation of traveling waves in nonlinear viscoelastic solids has also been 


;■ 0 , 0 , 0 , 


studied previously in the context of classical theory of viscoelasticity (see e.g. 

and references therein). The results of present work exhibit some similarities with those 

in the literature. The first common point is that the equations of motion admit kink-type 

traveling wave solutions. Also, in both cases, the effective width of the traveling wave 
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is proportional to the viscosity parameter and the wave profile becomes smoother as the 
viscosity parameter increases. However, our study differs from the articles within the 
context of classical viscoelasticity theory in the sense that the governing equation in our 
model (see (1.1)) is in terms of the stress and also the nonlinearity is on the inertia term. 

The structure of the paper is as follows. In Section[2]we introduce the one-dimensional 
strain-limiting viscoelasticity model as well as give a list of four nonlinear constitutive 
relations that has been suggested for elastic solids. In Section [3] we consider traveling 
wave solutions of the governing equations. In Section [ff] we solve the resulting differential 
equation for different nonlinear constitutive relations, and give analytical solutions where 
possible, or obtain numerical solutions. 


2. One-dimensional strain-limiting viscoelasticity 

Consider a one-dimensional, homogeneous, viscoelastic, infinite medium exhibiting 
small strains for large stresses. In the absence of external body forces, the equation of 
motion is given by 

PqUu = T x , (2.1) 

where po is the mass density of the medium, the scalar-valued function u(x,t) is the 
displacement, and T(x,t) is the Cauchy stress. Here and throughout this work the 
subscripts denote partial derivatives. In contrast to explicit constitutive relations of the 
classical theories of viscoelasticity, we shall employ an implicit constitutive relation 

e + ue t = g(T), (2.2) 


which gives the linearized strain e = u x and the strain rate e* as a nonlinear function of 
the stress T, with g( 0) = 0 and a nonnegative constant v. The model defined by (12.21) 
is the one dimensional form of When v = 0, it reduces to the one-dimensional 


version of the model introduced by Rajagopal in 




for elastic solids. 


For convenience, we now define the dimensionless quantities 


x = 


V 


t = 


t 


T = 


T 


u = 


v p 


V 


(2.3) 


where I is a characteristic length and p is a constant with the dimension of stress. 

Differentiating both sides of m with respect to x, substituting (12.21) into the resulting 
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equation and using dia , we obtain (HU), where we drop the overbar for notational 
convenience. The question that we shall discuss throughout the rest of this work is 
which of the possible forms of the nonlinear function g(T ) are relevant for the existence 
of traveling wave solutions of m- Following mainly the standard techniques used 
widely in the literature to find traveling wave solutions we obtain the solutions of ( 0 >, 
explicitly, implicitly or numerically, for various forms of g(T). 

We now discuss some strain-limiting models reported in the literature for elastic and 
viscoelastic solids. The following is a list of nonlinear constitutive relations g(T) which 
we adopt in this study. 

Model A : We first consider the one-dimensional version of the model proposed in 
an elastic setting by Kannan, Rajagopal and Saccomandi in Jsj], namely, 

5 (T) = /3T + a(l + !T 2 ) n T, (2.4) 


where a > 0, /? < 0, 7 > 0 and n are constants. Note that when n = 0 and/or 7 = 0, 
one recovers the standard constitutive equation for a linearized material. In Section [I] 
for the strain-limiting viscoelastic model defined by m-m we obtain traveling wave 
solutions explicitly if n = 1 and implicitly if n = — 1 / 2 . 

Model B : The second model is based on a simplified version of the nonlinear con¬ 
stitutive relation proposed by Rajagopal in Q: 


g(T) = 


T 


(1 + |Tj r ) 1 / r ’ 


(2.5) 


where r > 0 is a constant. This model was studied in elastic settings by many authors 
in different contexts (see e.g. m 3). Note that when /3 = 0, n = —1/2, a = 1 and 
7 = 2, Model A becomes equivalent to Model B with r = 2. In Section 0J when r = 2, 
traveling wave solutions corresponding to this model are obtained in closed form. 

Model C : This model is the one-dimensional form of the constitutive relation pro¬ 


posed by Rajagopal in 


BQ; 


g(T) = a 


1 — exp — 


PT 


jT_ 

l + 5\T\J\ ' 1 + |T| 


+ 


( 2 . 6 ) 


where a, /3, 7 and 6 are constants. Note that when /3 = 0 and a = 7 = 1 this model 
reduces to Model B with r = 1. In Section 01 we solve the nonlinear differential equation 
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corresponding to this model numerically and compute traveling wave solutions for a 
specific set of parameter values. 


Model D : This model is the one-dimensional form of a different model again intro¬ 
duced by Rajagopal in 0 0 ; 


g{T) = a 1 - 


1 


1 + 


+ 0 1 + 


1+<5|T| , 


1 


I+7T 2 


T, 


(2.7) 


where a, /3, 7 and S are constants. Note that when a = 0, with appropriate choice of the 
remaining parameters, we may derive Model A from this model. In Section [H traveling 
wave solutions corresponding to this model are also obtained numerically for a specific 
set of parameter values. 

Before going further, we would like to recall the remark made by Rajagopal in 0] 
about Models C and D. He says that both Model C and Model D have a drawback when 
the stress is compressive and sufficiently large. It is obvious from m and m that 
the assumption of small strain will be violated due to the initial terms in these equations 
when the stress is negative and sufficiently large. Furthermore, as it was also mentioned 
by Rajagopal in the same article, there are typographical errors in equations (3.12) and 
(3.13) of [141. where Models C and D were introduced. 

Figure Q] shows the variation of g{T) with T in a moderate stress regime for the above 
mentioned four nonlinear models with some specifically chosen parameter values. We 
note that, due to (12.211 . the vertical axis in Fig. |T] measures the sum of the linearized 
strain and the strain rate. Moreover, we observe that in the case of moderate stress 
levels the linearized strain may remain finite for the above models of strain-limiting 
viscoelasticity depending on the parameter values. 


3. Traveling wave solutions 

In this section we investigate traveling wave solutions of m- Traveling waves are 
solutions of the form 

r = T(0, z = x-ct, (3.1) 

where the wave propagation speed c is a constant to be determined below. Substitution 
of (13.11) into (11.11) reduces the third-order partial differential equation to the third order 
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Figure 1: Variation of the function g(T) with moderate T values for various nonlinear 
models of strain-limiting viscoelasticity. The specific sets of the parameter values are as 
follows: a = 0.5, /? = —0.01, 7 = 1 and n = —0.5 for Model A, r = 2 for Model B, 
a = 0.5, /3 = —0.01 and 7 = 5 = 1 for Model C and a = 0.5, /? = —0.01, 7 = 5 = 1 and 
n = 0.5 for Model D. 

ordinary differential equation in the variable £ given by 


T" —vc T'" = c 2 [g(T)]" 


(3.2) 


where the symbol ' stands for differentiation. For the rest of this study, we focus on 
traveling wave solutions of (13.21) that correspond to the heteroclinic connections between 
two constant states. Obviously, T(£) = constant is a trivial solution of (13.21) . so we 
assume that 


lim T(t) = T~, lim T(£) = T+ 


(3.3) 



with T ^ 7 ^ Tj+, where T ^ and T+ are constants to be specified later. Our main problem 
is to find restrictions on the nonlinear function g{T), which guarantees the existence of 
such a traveling wave solution, and is to discuss, from this point of view, the constitutive 
functions suggested in the literature. 

We now integrate (E2D once and then use the boundary conditions T'(£),T"(£) —> 0 
as £ —> ± 00 , to eliminate the arbitrary integration constant. A further integration of the 
resulting equation yields 


T-vcT' = c 2 g(T) + A, 
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(3.4) 









where A is an arbitrary integration constant. Boundary conditions m then give 


A=- 2 {T-+T+- c 2 [g(T~) + g(T+)}}, 


(3.5) 


and 


c 2 = 


T~ - T+ 

OO OO 


(rr-\_ (t+Y ( 3 - 6 ) 

9 {Too) 9 \Too) 

Thus the squared wave speed is obtained in terms of the two known states at infinity. 
Using (13.51) to eliminate A in (13.41) we get the differential equation 


T' = /(T) 


(3.7) 


where 


/(T)=— T- 


VC 


T - +T + 


g(T)~ 


g(T~)+g(T+) 


Of course, two obvious equilibrium points of (13.71) are T = and T = T+, that is, 
/(T-) = /(T+) = 0. Integrating (13.71) we get the implicit solution in the form 


£ - £o = 


'To 


ds 

7(a)’ 


(3.8) 


where £o is a constant and T(£o) — To- 

We conclude this section with a description of a prototype problem on which we will 
discuss the consequences of various forms of the nonlinear function g{T) in the next 
section. Recall that a heteroclinic traveling wave propagates from one constant state to 
the other if c 2 > 0. Due to (13.61) this implies that a traveling wave solution of (13.2D - (13.31) 
exists in one of the following two cases: 


Case (i) > T+ and giT^) > g{T +), (3.9) 


or 

Case (ii) T~ < T+ and g(T~) < g(T+). (3.10) 

For the remainder of this paper, without loss of generality, we restrict our attention to 
the first case for tractability reasons. Furthermore, we assume that the two constant 
equilibrium states are a normalized state of the stress and the zero reference state of the 
stress; that is, we take 

T~ = 1 and T+ = 0. 
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(3.11) 








One should recall that the stress is dimensionless. Also, even though we restrict our 
attention to the case (13.111) we should be aware that the nonlinearity amplifies the values 
of g(T ) when |T| > 1 and reduces them when |T| < 1, and, depending on which range of 
T we are working in, the traveling wave profile for g{T) is affected correspondingly. 

We note that (13.91) . (13.111) and g(0) = 0 imply g( 1) > 0. This condition is automat¬ 
ically satisfied by Model B (recall that r > 0), but it imposes the following restrictions 
on the parameters of Models A, C and D; namely, 


g( 1) = a 


5(1) = P + « (l + 70 >0 

P 


1 — exp I — 


and 


5(1) = 


1 + 6 
2 + 7 


2 + 6 V! + 7 

respectively. Plugging (13.111) into (13.611 gives 


2»°. 


> 0 , 


c 2 = 1/ 5 (1) 


(3.12) 


where we have used 5 ( 0 ) = 0. With the use of (13.111) . the differential equation (13.71) 
becomes 

T '=-^m mT - 9(T)] ' (3 - 13> 

which is studied for various forms of g(T) in the next section. We also note that solutions 
of (13.71) and (13.131) are translational invariant. That is, if T(£) is a solution of (13.71) or 
(13.131) . then so is T(£ + p) for any fixed constant p. Consequently, noting that T(0) 
can take any number in the range of values for T, we fix the traveling wave solution by 
assuming that 

T(0) = 1/2. (3.14) 

Two equilibrium points of (13.131) are clearly T = 1 and T = 0 (since g(0) = 0). 
Equation (13.131) may have additional equilibrium points depending on the form of g(T). 
Assume that (13.131) has an equilibrium point T* for which g(l)T* = g(T*). The lin¬ 
earization of (13.131) at this point possesses one real eigenvalue 

g{l)-g'{T*) 

vcg(l) 

which shows that T* is an unstable equilibrium for g( 1) ^ g'(T*), and a stable equilibrium 
for g(l) = g'(T*). 
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4. Applications to some nonlinear models 


This section discusses in detail both quadratic and cubic models of strain-limiting 
viscoelastic solids and the nonlinear models presented in Section 2, within the context of 
Section [3] 

We first remark that there is no heteroclinic traveling wave solution when we consider 
an elastic solid for which v = 0. This can be easily seen from EU or (13.131) by neglecting 
the derivative term (i.e. the dissipation term). Then, the only solution of the resulting 
algebraic equation is a constant solution but the boundary conditions at infinity require 
two different constants, giving a contradiction. 

A similar conclusion is also valid for the linear viscoelastic model for which we have 
g(T) = g'(0)T with g'( 0) ^ 0. In such a case (13.71) (or (13.131) ') reduces to T' — 0 which 
implies T is a constant. Following the same line of reasoning we find that there is no 
heteroclinic traveling wave solution for the strain-limiting linear viscoelastic model. 

In the remaining part of this section we focus on six particular forms of g[T)\ the 
quadratic and cubic models, and the nonlinear models described in Section (2j namely 
Models A, B, C and D. Figure [2] shows the variation of g(T) with T for linear, quadratic 
and cubic models in a moderate stress regime. We observe from Fig. [2] that, depending 
on the chosen parameter values, the quadratic and cubic models exhibit qualitatively dif¬ 
ferent responses and they may give rise to negative or large positive values of g(T) with 
increasing values of the stress. Obviously, the case where a positive (tensile) stress gives 
rise to a negative (compressive) strain is physically unacceptable in one-dimensional elas¬ 
tic or viscoelastic medium. Additionally, for large and positive values of g(T) the small 
strain assumption of strain-limiting viscoelastic solid is violated. Therefore, we conclude 
that, in general, the quadratic and cubic models may result in either physically unac¬ 
ceptable strain values or strain levels that are not consistent with the linearized strain 
assumption of strain-limiting theories. However, since they are the simplest representa¬ 
tives of the nonlinear models, for completeness we begin our discussion by considering 
general quadratic and cubic constitutive relations. 
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T 

Figure 2: Variation of the function g(T ) with moderate T values for the linear, quadratic 
and cubic models of strain-limiting viscoelasticity. The specific sets of the parameter 
values are as follows: g'(0) = 1 for the linear model, g'(0) = 1 and g"{0 ) = —0.6 for the 
quadratic model, and g'(0) = 1, g"(0) = —1 and g"'( 0) = 0.5 for the cubic model. 


4-.1. Quadratic case 

For the quadratic case, we assume that the function g(T) is of the form 

g(T)=g'(0)T+±g"(0)T 2 . (4.1) 


We first consider the traveling wave problem with the boundary conditions (13.31) and 
then specify them to be as in (13.111) . Substitution of (14.11) into (13.61) gives 

-i 


c 2 = 


ff'(0) + ^"(0)(T-+T+) 


The restriction c 2 > 0 requires that one of the following two cases must hold: 


fl"(0) > — 


-V( o) 


Too + Too 


and 


T~ + T + > 0 

^ OO ' ^ OO ^ 


or 

9"( 0)<-^ and T-+T+<0. 

With the use of (ED in ED, the differential equation we need to solve becomes the 
Riccati differential equation 

T' = B 2 T 2 + a\T + dp, 
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(4.2) 











where ag, ai and 02 are constants defined by 


a 2 = 


cg"{0 ) 

2u ’ 


_ 1 -cV(0) ^ 1 (rr + , n n\ 

a l — - — -) a 0 — —^ (^oo + ^oo) (1 ~ 


with 


1 „///^ (r~) 2 +(T+) 2 


0 = 


5 '( 0 ) + ± 5 "( 0 )^ 


+T~ 


g'(0) + i ff "(0)(T^ + T+) 

We observe that, when = 0 or T+ = 0, the coefficient ag vanishes and (14.211 reduces 
to the Bernoulli differential equation. This makes it possible to find explicit solutions. 




(a) Stress (b) Strain 

Figure 3: Variation of (a) T(£), and (b) g(T (£)) of the quadratically nonlinear model 
with £ for three different values of v (where g'{ 0 ) = 1 and g”{0) = — 0 . 6 ). 


We now turn to the case of (13.111) in which the wave propagation speed and the 
constants ag, a± and a 2 reduce to 


c 2 = 


</( 0) + \g"{ 0) 


> 0 , 


and 


g"{ o)c 

a 2 — ---, ai — — 02 , ag — 0 , 

2v 


respectively. Consequently, US) takes the form 

T' = a 2 T{l — T), 


(4.3) 


(4.4) 
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which admits the only two equilibrium solutions T = 0 and T = 1. Using m we find 
that, under the condition (13.141) . the explicit solution of (14.41) is found as 

T(0 = (l + exp(a 2 £)) ■ ( 4 -5) 

The important point to note here is that 631) satisfies the conditions defined by m and 
(13.111) if 02 > 0. Combining this with (14.31) implies that the traveling wave solution exists 
if g"{ 0) < 0 and c > 0 or if g"{ 0) > 0 and c < 0. In other words, the heteroclinic wave 
solution is a right-going traveling wave if g"(0) < 0 and a left-going wave if g"(0) > 0. 
On the other hand, there is no heteroclinic traveling wave solution for the quadratic 
model if 02 < 0, that is, if g"( 0) and c have the same sign. Figure [3] shows the variation 
of the analytical solution given in (14.51) with £ for three different values (corresponding 
to small, moderate and large values) of the viscosity parameter zz, as well as the profile 
for the sum of the linearized strain and the strain rate. We observe from Fig. [3] that, 
as it is expected, the traveling wave profiles become smoother as the viscosity increases. 
We also deduce from Fig. [3b] that the profile for g(T) is strongly distorted, in fact its 
values are reduced, due to the nonlinear dependence. Of course, this distortion can be 
intensified by choosing the values of the parameter g'{0 ) appearing in the constitutive 
relation properly. Furthermore, by choosing this parameter sufficiently small one could 
stay in the regime where the linearized strain assumption is valid. 

We close this part by examining the possibility of a shock wave (a traveling disconti¬ 
nuity). Differentiating the explicit solution T(£) given by (14.51) . we get 

no = ( 4 . 6 ) 

(l + exp(a 2 £)) 

The effective width of the traveling wave is defined as d = {T~— T+)/max |T'(£)|. Using 
(13.111) we conclude from (14.51) that the effective width of the heteroclinic traveling wave 
is d = 4/|a21 for the quadratic model. From (14.31) . it follows that d = /\g" (0)c\. Since 
the width d is proportional to the viscosity parameter zz, it is natural to expect that the 
wave profile becomes smoother as v increases. Furthermore, since the denominator of 
T'(£) in (14.61) is never zero, we conclude that a shock wave does not form. 
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4-2. Cubic case 


In this case we assume that 


g(T) = g'(0)T + ^"(0)T 2 + ig"'(0)T 3 . 


(4.7) 


Substitution of m into (13.71) yields 
1 


T' = 


T- 


T+ + T~ 

oo 1 oo 


g'(0)T + -g"(0)T 2 + —g'"{0)T 3 — 


g(T+)+g(T~) 


(4.8) 


We again consider the case defined by (13.111) . So it follows from (13.121) that 


c 2 = 


s'(0) + ^"(0) + ^' , '(0) 


1 -1 


> 0 . 


Using (13.111) in (14.81) (or using (14.71) in (13.131) ') we obtain the differential equation 

T’ = aT{l-T)(T+b), (4.9) 


where the constants a and b are given by 


a = — - 


cg"'{ 0) 


and 


6=1 + 3 


9"( 0) 


(4.10) 


Qv ” * ' g'"(0 )' 

It is worth pointing out that CD2) admits three equilibrium solutions; T = 0, T =1 and 
T = —b. Solving the differential equation (14.91) with (13.141) gives the closed-form solution 

rpl+b 1 


(1 -T) b (T + b) 1 + 2 b 


exp(6(l + b)a£). 


Figure U presents the variation of T(£) and g{T(£)) with ^ for three different values of 
<?"(0), namely, g"(0) = 0,0.25,0.75, when v = 0.5, g(0) = 1 and g"'{ 0) = 0.5. Note the 
different scales for the vertical axes of Fig. [4a] and Fig. [4b] We observe that the profiles 
for the stress and the strain become smoother as g"{0 ) (or equivalently b) increases. 
Similar to the quadratic case the profile for g(T) is distorted due to the nonlinearity. 
Also, as it is expected from the behaviour of the cubic nonlinearity in Fig. [2] the values 
of T are amplified in Fig. I4bl Furthermore, we note that, when g"( 0) = 0 (or equivalently 
b = 1), it is possible to obtain an explicit solution from the implicit one as 


no 


exp (aQ 

(3 + exp(2a^)) 1 / 2 


(4.11) 
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(a) Stress (b) Strain 

Figure 4: Variation of (a) T(£), and (b) of the cubically nonlinear model with £ 

for different values of </'(0) (v = 0.5, g'( 0) = 1 and g"'(0) = 0.5). The curves on the graph 
correspond to three different values: g"( 0) = 0, 0.25, 0.75 (or equivalently 6=1,2.5,5.5). 

The crucial fact is that the conditions given by (El impose a restriction on the constant 
a, which is the condition a < 0. This implies that, for the special cubic model, the 
traveling wave solution exists if a < 0, or equivalently, if g"'( 0) and c have the same sign. 
The variation of the analytical solution given in (14.111) with £ for three different values 
of the viscosity parameter v produces a figure, which is very similar to Fig. 3 and shows 
that the same conclusions are also valid for the present case, and therefore we do not 
reproduce it here. 

4-.3. Case of Model A 

We now take El to define the constitutive relation for our strain-limiting viscoelastic 
solid through (12.21) . Of course, if n = 0, El gives the linear model of strain-limiting 
viscoelasticity, for which we have already mentioned that there is no heteroclinic traveling 
wave. In general, depending on the values of the parameters appearing in (12.41) the 
function g(T) exhibits very different patterns of behaviour. As stated before, when 
/? = 0, n = —1/2, a = l and 7 = 2, Model A becomes equivalent to Model B with 
r = 2 and the implicit solution corresponding to this special case is given in the next 
subsection. In this subsection we restrict our attention to the case n = 1, which allows 
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us to find an explicit solution for the corresponding differential equation. Substituting 
(12.dl) with n = 1 into (13.131) yields 


T' = kT(1 — T 2 ), 


(4.12) 


where 


ay 


[a(l+ 7 ) + P\vc ^ ^ 

It is clear that (14.121) is a special case of (14.91) . with 6 = 1 and a = n. Therefore, if we 
replace a in (14.111) by k we get the explicit solution corresponding to (12.41) with n = 1. 
Additionally, we conclude that the traveling wave solution exists if k < 0, or equivalently 
if c[a(l + 7 ) + 0\ <0 (recall that v > 0, a > 0 ,7 > 0). Since the traveling wave solution 
is a special case of (14.111) , we can draw the same conclusions by saying that the traveling 
wave profiles become smoother as the viscosity increases as well as the wave profile for 
g(T) is distorted due to the nonlinearity. 


4-4- Case of Model B 




% 5 

(a) Stress (b) Strain 

Figure 5: Variation of (a) T(£), and (b) g(T(£)) of Model B with f for three different 
values of v (r = 2 ). 

Here we take (12.51) as g(T) in (12.21) . Note that combining (12.51) with (13.121) gives 
c 2 = 2 1 / r . If we substitute (12.51) into (13.131) . we get the differential equation that we need 
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to solve as 


T 


2 1 / r 


T' = — 1 - 

vc V (1 + |T| r ) 1 / r 


In the case of r = 2, we find an analytical solution of (14.141) . Using 
obtain the solution implicitly as 

H(T) = H( l/2)exp(£/i/c), 

where the function H(s) is defined as 

( I -* 2 ) 2 


(4.14) 
and (13.141) we 

(4.15) 


H(s) = 


(1 + s 2 ) 1 / 2 + 1 


^s[3 + s 2 + 2 3 / 2 (1 + s 2 )] ) \ s 

Note that some basic properties of H(s) are as follows: 


2 1 / 2 


-ff(Tl) = 0, 


H{ 1/2) > 0, 


and 


H(s) —>00 as s —y 0 + . 


By combining these properties and (14.15|) . we deduce the following two sets of results, 
depending on the sign of c. If c > 0, we get 

T ->■ 0+ as £ —>• +00 and T —> =f 1 as £ —► — 00 . 

Similarly, if c < 0, we get 

T y -pi as £ —y -Poo and T —y 0^" as £ —y 00 . 


We restrict our attention to the case c > 0 since the corresponding conditions are com¬ 
patible with (13.31) and (13.111) . 

Figure [5] shows the graph of the implicit solution for a specific set of parameter values 
and for three different values of the viscosity parameter u, and also the profile of the 
strain. From Fig. [5j we can see that the traveling wave profile becomes smoother as the 
viscosity increases, and the profile for g(T) is distorted due to the nonlinearity as in the 
previous models. 
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(a) Stress 


(b) Strain 


Figure 6: Variation of (a) T(£), and (b) g(T(£)) of Model C with £ for three different 
values of v (a = 0.5, j3 = 0.01 and 7 = 5 = 1). 


4-5. Case of Model C 

We obtain the differential equation we need to solve by substituting (12.61) into (13.131) . 
However, since the resulting equation is highly nonlinear and an analytical solution is 
not available, we focus on the numerical solution instead. To this end we use MATLAB 
function ode45 to solve the differential equation, which is the standard solver of MATLAB 
for ordinary differential equations. Omitting the details of the numerical calculations, we 
show in Fig. [6]the numerical solutions for three different values of the viscosity parameter 
v. 


4-6. Case of Model D 

When we substitute m into (13.131) we again get a highly nonlinear differential 
equation for which an analytical solution is not possible to find. For this reason, it is 
convenient to solve it numerically using MATLAB function ode45 just as above. In Fig. 
0we plot the numerical solutions for three different values of the viscosity parameter v. 

Both Fig. [6] and Fig. Q] clearly show that Models C and D have kink-type traveling 
wave solutions and that the wave profiles obtained numerically for the stress are in 
qualitatively good agreement with those derived from the analytical solutions belonging 
to the previous models. We note that the wave profiles for the strain are significantly 


18 














different from those of the previous models. The remarks made for those models regarding 
the smoothness and distortion of the profiles are also valid in both cases. We also observe 
that the reduction of the values of g(T ) in Model D is significantly stronger than that of 
Model C. 



% 



% 


(a) Stress 


(b) Strain 


Figure 7: Variation of (a) T(£), and (b) g(T{£)) of Model D with £ for three different 
values of v (a = 0.5, j3 = 0.01, 7 = <5 = 1 and n = 0.5) . 
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